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Abstract 

We argue that the massive Yang-Mills model of Kunimasa & Goto, Slavnov, and 
Cornwall, in which massive gauge vector bosons are introduced in a gauge-invariant 
way without resorting to the Higgs mechanism, may be useful for studying diffractive 
scattering of strongly interacting particles. With this motivation, we perform in this 
model explicit calculations of S-matrix elements between quark states, at tree level, 
one loop, and two loops, and discuss issues of renormalisability and unitarity. In par¬ 
ticular, it is shown that the S-matrix element for quark scattering is renormalisable 
at one-loop order, and is only logarithmically non-renormalisable at two loops. The 
discrepancies in the ultraviolet regime between the one-loop predictions of this model 
and those of massless QCD are discussed in detail. In addition, some of the similari¬ 
ties and differences between the massive Yang-Mills model and theories with a Higgs 
mechanism are analysed at the level of the S'-matrix. Finally, we briefly discuss the 
high-energy behaviour of the leading order amplitude for quark-quark elastic scatter¬ 
ing in the diffractive region. The above analysis sets up the stage for carrying out 
a systematic computation of the higher order corrections to the two-gluon exchange 
model of the Pomeron using massive gluons inside quantum loops. 


1 Introduction 


The quantitative description of diffractive phenomena within the framework of QCD is a 
long-standing problem. In part, the difficulty arises because diffractive processes involve 
both hard and soft scales resulting in a complicated interplay between perturbative and 
non-perturbative effects. One way to tackle this problem is to attempt a description using a 
“dressed” version of the perturbative degrees of freedom, where the “dressing” is meant to 
mimic the role of non-perturbative effects. Following Low [OJ and Nussinov ||, Landshoff 
and Nachtmann (LN) || introduced a two-gluon exchange-model of diffractive scattering, 
where they assumed that the infrared behaviour of the (abelian) gluon propagator is mod¬ 
ified by non-perturbative effects. Their success in reproducing several of the features of 
Pomeron exchange suggests that such an attempt may not be totally futile, and makes the 
question of how to compute systematically higher order corrections within this model all 
the more interesting. [] 

In the LN picture of the Pomeron the need for modifying the gluon propagator arises 
as follows: The simplest Feynman diagram which can model the Pomeron (exchange carry¬ 
ing the quantum numbers of the vacuum) is a box-diagram where two off-shell gluons are 
exchanged between incoming and outgoing quarks, which scatter elastically. The pertur¬ 
bative calculation of the above process gives rise to an infrared singularity at t — 0, whose 
origin is the fact that the bare gluon propagator do(q 2 ) diverges at q 2 = 0. Specifically, the 
amplitude obtained from such a diagram assumes the form where /3q is 

given by /3 q ~ f dq 2 [a s d(q 2 )} 2 , and d(q 2 ) is the gluon propagator. The introduction of a 
“massive” gluon propagator is the simplest way to obtain a finite /3 q and the gluon mass is 
then fixed by data. 

It has been suggested long ago 0 that the non-perturbative dynamics of QCD lead 
to the generation of a dynamical gluon mass [] while the local gauge invariance of the 
theory remains intact. This gluon “mass” is not a directly measurable quantity, but must 
be related to other physical parameters such as the string tension, glueball masses, or the 
QCD vacuum energy 0, and furnishes, at least in principle, a regulator for all infrared 

* A different approach is provided by the BFKL formalism Q], which is the most serious attempt at 
a first principles QCD derivation of Pomeron exchange to date. However, the perturbative nature of the 
BFKL approach often makes it unsuitable for the analysis of diffractive scattering, where both soft and 
hard momentum scales are in general relevant. 

^Dynamically generated masses depend non-trivially on the momentum; in particular, they vanish for 
large momenta. This property is crucial for the renormalisability of the theory ||. 
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divergences of QCD. The above picture emerged from the study of a gauge-invariant set 
of Schwinger-Dyson equations In addition, lattice computations || reveal the onset 
of non-perturbative effects which can be modelled by means of effectively massive gluon 
propagators. Various independent held theoretical studies spanning almost two decades 
TOPI also corroborate some type of mass generation, although no consensus about the 
exact nature of the mass-generating mechanism has been reached thus far []. Interestingly 
enough, the effective gluon propagator derived in 0 describes successfully nucleon-nucleon 
scattering when inserted, in a rather heuristic way, into the two-gluon exchange model [T3| . 
Despite this phenomenological success however, it is not clear whether a dynamical gluon 
propagator may be used in calculations as if it were a tree-level propagator derived from 
Feynman rules. More importantly, it is not known how to systematically improve upon 
such a calculation, i.e. how to compute higher order corrections. 

Given this lack of a computational scheme originating from a “first principles” QCD 
treatment, we propose instead to resort to a field theory which is formally close to QCD 
and contains at the same time the feature which appears to be phenomenologically use¬ 
ful, namely a gluon mass. To that end we revisit a model introduced independently by 
Kunimasa & Goto [[L4|, Slavnov [|15|], and Cornwall JR|, which provides the extension to a 
non-Abelian context of the work of Stueckelberg [|17]]. This model accommodates massive 
vector bosons without compromising local gauge invariance and without introducing a Higgs 
sector. In what follows we will refer to it as the massive Yang-Mills (MYM) model. 

In the MYM model, a mass term is added directly to the Yang-Mills (YM) Lagrangian 
and gauge invariance is preserved with the help of auxiliary scalar fields. Unlike the usual 
Higgs mechanism |L2| however, there are no additional physical particles appearing in the 
spectrum (no Higgs bosons). The price one pays is that perturbative renormalisability is 
lost. In particular, the one-loop 5'-matrix element for gluon elastic scattering gg —> gg is 
known to be non-renormalisable ]llj; its renormalisability can be restored only with the 
introduction of Higgs boson in the spectrum [|l^^,^l|] . This fact renders the MYM non- 
renormalisable at one loop. However, the introduction of a Higgs boson is not necessary for 
the renormalisability of the one-loop 5'-matrix of the processes qq —> qq which is relevant 
for diffractive scattering. As we will see in detail, the first time this latter process receives 
(logarithmically) non-renormalisable contributions is at two loops. In addition, the model 
has been shown to be unitary, in the sense of the optical theorem, to all orders in pertur- 

*Of course, the introduction of a gluon mass at tree level through the usual Higgs mechanism |l^] is 
excluded, as it would introduce extra (unwanted) scalar particles in the physical spectrum. 
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bation theory |H|. Several formal properties of this model have been extensively studied 
in the literature cited above and are well understood. 


Our main phenomenological motivation for turning to the MYM model is to carry out 
the next-to-leading order corrections to the two-gluon exchange process for qq —> qq, in the 
context of a concrete field theory, where the effects originating from the presence of a gluon 
mass can be studied systematically. Clearly, before attempting such a complex calculation 
it is necessary to develop some familiarity with the predictions of the MYM model at 
leading and next-to-leading order. The purpose of the present work is to provide a detailed 
analysis of various field-theoretical issues which appear when one uses the MYM model 
for computing S-matrix elements | involving quarks as external states. In addition to the 
clarification of theoretical points, several of the results presented in this paper constitute 
useful ingredients of the full calculation. 

More specifically, we discuss the following points: 


• We analyse in detail how the MYM and QCD differ already at tree-level, and how this 
difference propagates to higher orders. In particular we show using both unitarity 
and analyticity arguments as well as explicit one-loop calculations how the tree-level 
discrepancy affects the one-loop beta function, i.e. it alters the high energy behaviour 
of the theory. 

• We verify explicitly in the context of a specific example that the S-matrix contains 
no unphysical poles. The cancellation of such poles, which is expected from formal 
considerations, provides a non-trivial consistency check of the model, and can serve 
as a guiding principle when carrying out lengthy calculations. 

• We demonstrate that at the one-loop level the scattering amplitude of interest is 
renormalisable, and that one can construct a gauge-invariant running coupling (effec¬ 
tive charge) just as in QCD. This leads to the definition of a gauge-invariant gluon 
propagator, generalising Cornwall’s construction for the standard QCD case. A de¬ 
tailed comparison of our result with the QCD one is performed. 

s By working directly with S-matrix elements one has the additional advantage of avoiding pathologies 
which affect individual, unphysical Green’s functions. In fact, because of several cancellations taking place 
at the level of S-matrix elements, the final answer often has better properties than those of the Green’s 
functions involved in the calculation. A typical example of this situation arises when using the unitary 
gauges for the electroweak model; in these gauges, Green’s functions are non-renormalisable, while S-matrix 
elements are |22|. 
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• We show that the non-renormalisable contribution arising at the two-loop level de¬ 
pends only logarithmically on the cutoff. This result is new, to the best of our 
knowledge; its derivation relies crucially on extensive cancellations which take place 
at the level of the S'-matrix after the judicious exploitation of the tree-level Ward 
identities of the MYM. 


The paper is organised as follows: In Section we briefly review the MYM formalism, and 
establish connections which will be useful for the calculations which will follow. In Section 
we analyse qq annihilation into two gluons at tree level within the MYM model, and 
compare with the result in standard YM. In Section f| we study the one-loop contributions 
to qq —> qq and show in detail how the MYM model gives rise to renormalisable and 
unitary S-matrix elements. In Section || we turn to the two-loop contribution to qq —> qq, 
and demonstrate the emergence of logarithmically divergent non-renormalisable S'-matrix 
elements. In Section |6] we investigate quantitatively the connection of the MYM to held 
theories where the gauge bosons acquire masses by means of the usual Higgs mechanism 
1 % . In particular, we show how the presence of a Higgs boson cancels the logarithmically 


non-renormalisable contributions found in the previous section. Throughout Sections || - 
[6] we use the pinch technique (PT) rearrangement of the S-matrix p3|fl in order to make 
several cancellations manifest. We hasten to emphasise, however, that the PT only serves 
as a convenient intermediate step, helping to expose the unitarity and renormalisation 
properties of the S'-matrix, but none of the final results reported here depends on the use 
of this method. In section [7] we take a first look at a possible phenomenological application 
of the MYM model, namely, quark-quark elastic scattering in the diffractive region. Finally, 
in Section § we summarise our results and discuss possible future applications. 


2 The Massive Yang-Mills Model 


In this section we first review briefly how local gauge-invariance and massive gauge bosons 
can be reconciled in the MYM. Next we show that the MYM is physically equivalent to 
a held theory where the gauge bosons have been endowed with a mass “naively”, i.e. by 
adding a mass term at tree-level without preserving gauge invariance. 


In order to introduce the MYM model |L5|,[[(|], let us start from the standard YM 
action for the SU(3 ) gauge group: 

1 


S Y M[A\ = --j d‘x Tv(r^) 


( 2 . 1 ) 
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where T^x) = d fi A u (x) - d u A fl (x) + ig [A^x), A„(x)\ and A^x) = A“(x)T a , with T a 
the SU(3) generators in the fundamental representation. For the purpose of the present 
discussion, matter fields can be ignored. Under a gauge transformation, parametrised by 
U(x), A M —> Au where 

A u Jx) = U(x) A^x) U~ 1 (x) - -U(x)d lx U~\x). (2.2) 

9 

The requirement of gauge invariance for the action forbids a naive mass term for the gluon. 
However, by introducing 577(3)-valued fields, V(x), one can define 


C^x) = --V(x)d li V~ 1 (x). (2.3) 

Under a gauge transformation one postulates that V —> V u = UV. As a consequence, 
Cfjx) has the same gauge transformation properties as the gauge held A^, i.e. 


CV = UCJJ~ l - -UdM~\ 

^ ^ g M 


(2.4) 


The quantity 


B li [A,V]{x)=A li {x)-C li {x) (2.5) 

thus transforms as = UB^U^ 1 under a simultaneous gauge transformation of the A M 
and V helds, so one can add to the Yang-Mills action the following gauge-invariant term: 


S m [A, V] - M 2 I d‘x Tr 


( 2 . 6 ) 


More explicitly, gauge invariance of the above quantity can be written as 


Sm[A u ,V u )=S m [AV}. (2.7) 

Recalling (|2.3|) and (l2~5|) , it is clear that Sm generates a mass term for the gluon held 
a kinetic term for the held V, and an interaction term between the A and V helds. 

Finally, we can write down the gauge-invariant action functional for the MYM theory: 


Smym[A, V ] = Sym[A] + Sm[A, V], 


( 2 . 8 ) 


We write now the path integral for such a theory. Gauge invariance of Smym implies that 
a gauge-hxing prescription is needed to quantise the theory. The Faddeev-Popov procedure 
can be carried out as in standard YM theory, leading to 

Z = j VV VA e iSMYM[A W A[A] 5(g[A]). (2.9) 
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Here the gauge-fixing condition is g[A] =0 and A [A] is the corresponding Faddeev-Popov 
determinant . In order to make the theory amenable to a perturbative treatment one could 
rewrite Sm[A, V] as a power series in the coupling constant g. This is obtained by writing 

V(x) = exp (ig 9 a (x)T a ), (2.10) 


and inserting the power expansion for V into ( |2.6| ). The resulting expression contains 
interaction vertices with an increasing number of scalar 9 fields and zero or one gauge field 
Ap. Then, using standard techniques, Feynman rules can be derived |jl5|. However, as long 
as one is only interested in gauge-invariant calculations, a considerable simplification of the 
Feynman rules can be achieved. To see this, let us consider the calculation of the vacuum 
expectation value of a generic gauge-invariant operator 0[A]: 

1 


(O) = — [ VV VA e iSMYM[A ' v1 A[A] 6{g[A])0[A]. 

£/ J 


( 2 , 11 ) 


We perform a change of integration variable in the VA integral, i.e. we rewrite it in terms 
of a new field A , defined through the following identity: 


Ar = VA'„ r'-W-'EA;! (2.12) 

In other words, A M and A^ are related by the gauge transformation generated by V. Thus, 
VA = VA. Also, gauge invariance implies that Sym[A] = Sym[A ], 0[A] = 0{A] and 
A [A] = A [A']. Strictly speaking, the last equality holds only if one neglects the issue of 
Gribov copies. This is correct for the purpose of a perturbative treatment. 

The crucial observation is that because of ( |2.3|) , ( |2.4|) and ( p.5|) one has 


BJA,V] = B m IA' v ,Vl = VA' lt v~\ 


(2.13) 


hence we can write 


S m [A, V] = M 2 J Ax Tr A^A G 

The path integral (j2 .1 1|) can then be rewritten as 


(2.14) 


<°> s * 


VV VA e 


(Sym[A\+M 2 jd*x TV A'„AV) a[a , } s(s[a , V]) ^ p 15) 


Notice that in the above expression all the dependence on the V fields is carried by the 5- 
function. The integration on the V fields yields a factor 1/A [A'], which cancels the Faddeev- 
Popov determinant arising from the gauge-fixing procedure. The final path integral can be 
written as 


(O) = ^ J VA e Ks™[A]+M* f d*x Tr A^) ( 2 . 16 ) 
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The above manipulations show that, as long as the operator of interest is gauge invariant 
in the usual massless QCD sense, the model defined by Q2.9|) is equivalent to the simpler 
massive vector theory defined by ( |2.16|) . The latter is obviously much easier to handle in 
perturbative calculations. 

It is important to emphasise that the models are not equivalent at the level of (gauge 
dependent) Green’s functions of the gluon field. In particular, let us compare the tree-level 
expressions for the gluon propagator in the two models. From ( [2.9D one obtains (in the 
Landau gauge) 


while (|2.16| ) yields 


D^ e (k) = 

1 


kf.kA 

( 2 . 17 ) 

k 2 -M 2 

k 2 ) ’ 

DlT(k) = 

1 

( 9nv — 

k, t kA 

( 2 . 18 ) 

k 2 — M 2 

M 2 ) ' 


The former expression corresponds to a gluon with two polarisation states, as in the massless 
case, whilst the latter has three polarisation states, as expected for a massive vector boson. 
Of course the number of degrees of freedom in the two models has to match. In fact, the 
third polarisation state of ( [2.18|) corresponds to the massless scalar field 9(x) which appears 
in the MYM model. 

We have seen that gauge invariance can be maintained in a theory of massive gluons 
without introducing additional particles into the spectrum. However, as we will discuss 
later, and as noted by others [|T8|, the resulting theory is no longer renormalisable. 


3 Tree-level analysis 


q(pi ) 


g{h), g,a 


q(p2) 


<? 


(a) 




(c) 


Figure 1: Diagrams contributing to Tff. 

In this section we will study in detail the tree-level cross-section for quark-antiquark 
annihilation into two massive gluons, i.e. qq —> gg, within the framework of the MYM 
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model. The reason is three-fold: First we want to gain some familiarity with the formalism, 
second we want to study the difference between the MYM and standard QCD at the level of 
physical amplitudes, and third, in conjunction with the results of the next section, we will 
check explicitly that the MYM model produces unitary ^-matrix elements. Throughout 
this section we use the methodology and notation first introduced in ||24|| . 

Let us consider the quantity A, 

A = \j d(p.s. 2 ) (qq\T\gg)(gg\T\qtf , (3.1) 

where 

I d(P.S ., 2 ) = I d% J d 4 k 2 5+(k\ - M 2 )8+(kl - M 2 )S^\q -h- k 2 ) ( 3 . 2 ) 

is the phase space integral for two particles with equal mass M in the final state, with 
5 + {k 2 — M 2 ) = 9(k°)5(k 2 — M 2 ). In (fTT|), the factor 1/2 is statistical, arising from the fact 
that the final on-shell gluons should be considered as identical particles in the total rate. 
A is the contribution to the imaginary part of the amplitude for qq —> qq which arises from 
a gluon loop. We first focus on the tree-level amplitude T = (qq\T\gg). Diagrammatically, 
the amplitude T consists of two distinct parts: t and w-channel graphs that contain an 
internal quark propagator, 7)//, as shown in Fig.|I](a,b) and an s-channel amplitude, T s //, 
as shown in Fig.jl](c). The subscripts “s” and “t” refer to the corresponding Mandelstam 
variables, i.e. s = ( p\ + p 2 ) 2 = (&i + k 2 ) 2 and t — {p\ — ki) 2 = ( p 2 — k 2 ) 2 . 

Let us first define the following quantities: 



y C (Pl,P2> = Jt'(P2)r C 7 p !l(pi) , 

rcjV.K.g) = gf ak D 0 (q)v‘(p 1 , P2 ), 

(3.3) 

where 




M q ) - q 2_ M 2 ■ 

(3.4) 

The amplitude is given by 

/-j-’O'b _ f-r ab , /-r a b 

(3.5) 

with 

Zf^nfr^q-k^-h), 

(3.6) 

where 

^ \fiu{q,-ki,-k 2 ) = (q + h) u gx» + {k 2 ^ki) x g^ ~ (q + k 2 )^g X u, 

(3.7) 
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is the usual three-gluon vertex and 


%£ = -ig 2 v(j>2)(T l >'y''—!- jr T<‘Y + r 1 >'- r i-y'T‘)«(p,). (3.8) 

^ v Pi- fa p 1- fa ' 

Notice that in (|3.6| ) only the part of the tree-level massive gluon propagator appears, 
since any longitudinal part vanishes due to current conservation when it hits the external 
on-shell quarks. The three-gluon vertex satisfies the fundamental Ward identity: 


,{q,-fa, -k 2 ) = [d 0 1 (k 2 ) - d 0 1 (q)}g Xu + [q x q„ - k 2X k 2v \ 

= [D~ 1 (k 2 ) - D~ 1 (q)]g Xv + [q x q u - k 2X k 2v \ (3.9) 


(and cyclic permutations) where d 0 1 (q) = q 2 . The form of the Ward identity in the massive 
theory is therefore identical to that of the massless theory. 


We then have 


A = 


2 

1 


T‘ b Q^(h) Q* x (k 2 ) v»d(p.s:- 


I/IV 


%t +r i 


/- j-ab f . /T-afef 
1 s a \ -r 1 t. 


tcr\ 


d(P.S . 2 ), 


(3.10) 


where 

Ut-L i.v 

cr{k) = -<r + — ( 3 . 11 ) 

is the polarisation tensor of the massive gluon. On shell, i.e. k 2 = M 2 , we have that 
k^Qfivik) = 0. This fact motivates the standard PT decomposition of the three-gluon 
vertex P5[ |: 

T Xfl faq,-fa,-k 2 ) = T^ u (q,-fa,-fa) + T^faq,-fa,-k 2 ) (3.12) 

where 


r A^(<h -fa, ~k 2 ) = (fa - fa) x g^ v + 2q v g XjJ - 2q fl g x „ 

/luid: k\, k 2 ) k\^g Xu k 2u g Xfl . 

The term vanishes when it hits the polarisation tensors, and fl3.10|) becomes 


(3.13) 


A = 


faf + %$\ cr(ki) Q- X (k 2 ) [rjf t + T^Up.s 2 ), 


(3.14) 


where 


ry F,ab 

1 -spuV 


sjp abj^F 
'^'p PPLV’ 


(3.15) 
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To evaluate further the expression on the RHS of (|3.14|) and establish its connection to 


massless QCD we proceed to determine the action of the longitudinal momenta coming 
from Q ,J - a {k x ) and Q vX {k 2 ) on %^ b and T t ab v . 


UUg-F,ab 

K l 2 Sfj, u ~ 

Pi - k 2 ) x k 2u - M 2 g Xu ]n ab - D^\q)1l ab , 

(3.16) 

l„vq-F,ab _ 

h ' J 2 I Sf lu — 

Kh - k 2 ) x h, + M 2 gx ,]n ab + D^\q)n ab , 

(3.17) 

UHq-ab _ 

A 1 V 

Do\q)K\ 

(3.18) 

ivryab _ 

~ 

p 

1 0 

Cl 

1 

(3.19) 


The terms proportional to D 0 l (q) cancel when forming the sum kf[T s ^ff b + T t ab u \, giving- 
rise to 


+ = [{k^k 2 ) x k 2v -M 2 gl]R ab , 

k"[TFf + Ttf u \ = [(k 1 -k 2 ) x k lll + M 2 g x ]n ab . (3.20) 

Such a cancellation is instrumental for the good high-energy behaviour of the resulting 
amplitudes. Using the longitudinal momenta inside the polarisation tensors to trigger the 
identities listed above, we can decompose A into three parts: 

A = Ai + Ao + As (3.21) 

where 


Ai = 

L 

f [T, F Tp - R p { 7 -(k, - hr(k, - k 2 y + 2M 2 <T} Rl\d(P.Sp„ 

(3.22) 

a 2 = 


r (%Tp + Tfp)d(P.S. 2 ), 

(3.23) 

A3 = 

1 

2 ./ 

! T t pd(p.s:‘). 

(3.24) 


Ai contains the purely propagator-like (self-energy) contributions, A 2 contains the vertex¬ 
like contributions and *4. 3 contains the box-like contributions. We see that all terms pro¬ 
portional to M~ 2 or M~ 4 have disappeared. Therefore, at this point, it is clear that at 
the one-loop level the MYM model gives rise to a renormalisable S-matrix for qq —> qq, 
provided that we assume unitarity and analyticity (i.e. dispersion relations). In the next 
section we shall check this conclusion by an explicit one-loop calculation. 

We now focus on the propagator-like part, A\ . Current conservation allows us to 
make the replacement 

- 8+ 4 (h - k 2 ) p (k 2 - k 2 ) x . (3.25) 
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Then ( |3.22|) becomes 

M = g 2 c A D 2 „(q) v;| J [(V - M 2 )g <“■ + l(h- k 2 f(h - k 2 )‘ 


d(P.S?)\vy, (3.26) 


where Ca is the Casimir eigenvalue in the adjoint representation. The final step is to use 
the following results for the phase space integrals: 

I d(P.S , 2 ) = A. %»)0(« 2 - 4M 2 ) A(<j 2 ), 

Id(P.S. 2 )(k 2 - k 2 ) F (k 2 - **)„ = - 4A7 2 ) 9 2 A 3 ( 9 2 ) 9j ,„ , (3.27) 


where 


We obtain 


A(g 2 ) = Wl 


AM 2 


(3.28) 


A = B 2 ( g )^{^cVA(? 2 )lf+ 


= Dl( q )vA^ c A q 2 A(, 2 ) 


/II 


1M 2 


V 3 24 J ^ 2 g 2 

29 


<r\v; 


(3.29) 


with a s = g 2 /(An). The reason why we write the coefficient A as the deviation from -A on 
the second line of Eq. ( |3.2S|) will become clear in what follows. 


It is instructive to repeat the same calculation for the case of massless QCD, in 


order to examine the physical difference between the two theories at tree-level |24|. The 
crucial modification, in the case of QCD, is that in ( 13.101 ) the polarisation tensors Q 
corresponding to the massive gluons, are replaced by the polarisation tensors P^ iu (k,r]), 
given by 


if) g^ u + 


VtdK + VJ:,, _ k^k v 


(3.30) 


gk '' (jjk) 2 ’ 

which are appropriate for massless spin-1 gauge bosons. As before we have that, for massless 
on-shell gluons, k^P^ = 0. All other expressions can be obtained directly from the MYM 
expressions simply by setting M 2 = 0. In particular, both the derivation and the final form 
of the Ward identities of ( |3.20|) are identical |[24l , [26| 

The QCD expression corresponding to (|3.26| ) is given by | |24] 


AT D = .9 2 c -a d 2 (q) Vf {/ [Aq 2 g^ + (h - - k 2 y]d(P.S 2 ) }v„ c (3.31) 
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and, after carrying out the phase space integration for the two final (massless) gluons, we 
obtain the QCD analogue of ( |3.29| ): 


Ai° D =«(?)■ 


T (t) c ^ F’ 


(3.32) 


Notice that the factor -d- in (|3.32|) is the characteristic coefficient of the one-loop (3 function 
of quarkless QCD. 

Obviously, if we set M 2 = 0 in ( |3.26| ) and ( [j.29| ) we do not recover the massless QCD 
result, i.e. A\(M 2 = 0) 7 ^ A® CD . In that limit the two answers differ by the amount 
dp this descrepancy heralds the difference in the leading logarithmic behaviour of the two 
theories, which we will establish in the next section. On the other hand, it is clear that 
Ai(M 2 = 0) = Af CD for i = 2,3. Evidently, even though the two theories satisfy the 
same type of tree-level Ward identities, the fact that we have to use different polarisation 
tensors for massive and massless gluons gives rise to different 5 -matrix elements, and this 
difference persists even in the limit M —► 0. As explained by Slavnov Jl^], the physical 
reason why the limit M —> 0 of MYM does not recover massless Yang-Mills is that one 
cannot continuously go from three polarisation states to two. It is interesting to notice 
that after the PT rearrangement the discrepancy between the two theories as M —> 0 has 
been isolated in the universal, process-independent, propagator-like piece, A\. 


4 One-loop analysis 

In this section we turn to the issue of unitarity and renormalisability at one-loop. To begin 
with, we show using a one-loop calculation that all unphysical poles introduced by the 
gauge-fixing choice cancel in an S'-matrix element. This cancellation is a necessary condition 
for proving the unitarity of the resulting expressions; indeed, if expressions containing mixed 
poles had survived, they would give rise to unphysical thresholds. Next, by comparing the 
results of this section with those of the previous one, we will be able to establish explicitly 
the validity of the generalised optical theorem to lowest order and hence have an explicit 
demonstration of unitarity at one loop. Finally, we show that the resulting expressions 
can be made finite by the usual mass and wave-function renormalisation. Throughout this 
section we employ the PT, which makes cancellations particularly easy to track down. 

We study the one-loop amplitude, AA = {qq\T\qq), for the process qq —> qq, using the 
Feynman rules derived by Slavnov Jl5| : the massive gluon propagator in the Landau gauge 
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0 is given by, 


D^(k) 


1 ( k lt k u \ 

k 2 - M 2 \ g ' w k 2 ) ’ 

D 0 (k) t[iv , 


(4.1) 


the ghosts are massless and only appear inside closed loops (with a statistical factor —1/2), 
and the three- and four-gluon vertices are those known from massless QCD. Note that we do 
not include quark loops since they are trivially related to the equivalent QCD diagrams and, 
as such, need not be considered when investigating the new features of the MYM model. 
We will show explicitly that all unphysical poles (i.e. massless poles in the Landau gauge) 
induced by the longitudinal part of D^ and by the massless ghosts vanish in the one-loop 
amplitude AL Moreover, all contributions containing unphysical poles are propagator-like, 
in the sense defined by the PT re-arrangement of the amplitude [^3|,U . 

First we define 


I(q,k) = 
J{k) ee 


{k 2 - M 2 )[{k + q ) 2 - M 2 } ’ 

1 


k 2 -M 2 ’ 

and the auxiliary expressions containing mixed poles: 

1 


I 0 (q,k) = 

h(q,k ) = 
h (q,k) = 
Ji(k) = 


k 2 (k + q ) 2 


k 2 (k 2 — M 2 )[(k + q ) 2 — M 2 ] ’ 

1 

k 2 (k + q) 2 {k 2 - M 2 )[(k + q ) 2 - M 2 } 
1 


(4.2) 


(4.3) 


k 2 (k 2 — M 2 ) ’ 

which appear in intermediate steps but vanish in the final answer. In addition, we define 


U^(k) = D 0 1 (k)g[ lv - k^k v 


(4.4) 


We consider the diagrams of Fig.2 individually. The expressions for all non- 
propagator-like contributions are the same as the corresponding contributions of the regular 

^Slavnov’s choice of the Landau gauge was motivated by the fact that it leads to a reduction in the 
number of interaction vertices. Of course, for computations of S-matrix elements any other choice of the 
gauge fixing parameter will lead to the same final answer 
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(d) 



Figure 2: One loop diagrams. 

QCD graphs in the Feynman-’t Hooft gauge, with the only difference that the internal gluon 
propagators are Dq rather than d,Q. These results emerge at the end of a gauge-invariant 
calculation and are not linked to any particular gauge choice. Consequently we turn our 
attention to the propagator-like contributions. 

For each diagram of Fig.2 we write the associated amplitude as a sum of propagator¬ 
like (P) and non-propagator-likc (. NP) pieces: 

M® = Mf+M { $p, (4.5) 

where i labels which diagram is being considered. Of course -M^p = M^p — Mnp — 0. 
For the propagator-like piece it is convenient to write 

r r\ n h 

Mf = g 2 c A Dl(q) V?VZ J n® (4.6) 

For graph (a) one finds that 

i$> = ng„ + ng„ (4.7) 
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where 


(4.8) 


n oi = V„,(q,k,-k - q)T^(q,k,-k - q) 


and 


n (a) 

n i,^ 


= ~h(q,k) 


:h(q, k) 


Dq 2 (q)gfj,u - 2D^(q)U-\k + q) + U~Hk + q)U p ~\k + q) 


D^ 2 (q)+2M 2 D^{q)+M 


k/jK. 


(4.9) 


nf, a) is the part of (a) which arises due to the g px g aT part of t pX t ar and contains only the 
physical (massive) poles. The unwanted mixed poles live in n^*. We note that the term in 
parenthesis accompanying the J 2 factor of (|4.9|) is equal to q A . We choose not to simplify 
this expression since the terms proportional to inverse powers of Do are going to cancel 
against similar contributions from vertex and box graphs; retaining them explicitly will 
make the mechanism of the cancellation more transparent. 

We next turn to the vertex graph of Fig.2(b) (and its mirror image). We write 



(4.10) 


The contribution which arises from the term in the product of the two gluon 

polarisation tensors is equal to the usual QCD vertex graph in the Feynman-’t Hooft gauge 
with massive, instead of massless, internal gluon propagators. This term can still give a 
propagator-like contribution due to the pinching of the fermion propagator triggered by 
the three-gluon vertex [E3]. This contribution is 


<^ = 2 1{q,k) Do(q)- 1 ^ 


(4.11) 


The n[ b) term contains the remaining parts of the polarisation tensor product and the 
pinching of the quark propagator is triggered by the momenta therein: 


Wl = 2D^(q)I 1 (q,k) 


~D 0 ( q)h{q,k) 


D 0 \q)g^ ~ U \k + q) 


Do\q) + M 2 


k t ,k v . 


(4.12) 


For the box graph of Fig.2(c) (along with the crossed box): 

n= ~Do 2 (q) h(q, k ) g^ + ^Dg 2 (q ) I 2 (q, k) k^k u , (4.13) 
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and for the remaining graphs: 


n (d) 

fLV 

= D 0 1 (q) Ji(k)g^, 

(4.14) 

n (e) 

fLV 

= k 2 Ji(k) t lxu (k) , 

(4.15) 

n (/) 

fLV 

1 

= -2 J o (<?,&) Mr- 

(4.16) 


Notice that, at this point, all terms containing massless propagators are multiplied 
by inverse powers of D 0 (q). If we now add these contributions, all terms proportional to 
inverse powers of Do(q), and therefore all terms containing massless poles, cancel against 
each other. In order for this cancellation to go through it is crucial that the ghost diagram 
has a statistical factor of (—1/2), rather than the —1 of massless QCD. It is also interesting 
to observe that the aforementioned cancellations take place algebraically before any of the 
integrations over the virtual momenta k are carried out. In particular, we have not resorted 
to the use of dimensional regularisation results such as / = 0 or / = 0. 

Our final result for the propagator-like part of the one-loop amplitude is thus 

n /w = n^ + I(q,k)(2DQ 1 (q)g fU , - (4.17) 

Notice that the last term in the above equation could be interpreted as a contribution from 
massive ghosts. This term has emerged naturally from our calculation, even though we 
started out with massless ghosts. 

As we have already mentioned, in the limit M — » 0, the box-like and vertex-like 
parts of the MYM S'-matrix element exactly reproduce their massless QCD counterparts. 
The propagator-like piece is, as expected from the work of the previous section, different. 
Defining the effective gluon self-energy, II, via 

M P = D 2 0 (q)V c »n(q 2 )V„ c (4.18) 


a straightforward calculation yields 





(4.19) 


where e = 4 — n. Setting M 2 = 0, and using / ^ = 0, the above expression reduces to 


n(<? 2 ) 


M =0 


2 2 
= g c A q 



(4.20) 
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The corresponding gauge-invariant effective self-energy for massless QCD, first given 


in |25||, reads: 


fi QCD (g 2 ) = 9 2 c A q- 2 


11 


+ 


d n k 


■h (q,k). 


(4.21) 


3 6/ J i(27r) r 

Notice that the above result can be obtained from (|4.17|) by taking the M 2 —> 0 limit and, 
at the same time, changing by hand the coefficient in front of the ghost term from (-1/2) 
to (-1). 

To establish contact with the previous section, we need to compute the imaginary 
part of M.p. Using 


3m 


d n k 

?(27t) 7 


I(q,k) 


3m <; 


| J dxln[M 2 — q 2 x(l — x)]j 
A (q 2 ) 

J d(P.S 2 ), (4.22) 


167T 2 >- JO 

6(q 2 - 4M 2 ^ 


167T 


it is straightforward to check that unitarity holds, i.e. 


23mA! p — A\. 


(4.23) 


Similarly, one can demonstrate the unitarity of the vertex- and box-like contributions. 

We next proceed to renormalise the expression for II(g 2 ); we carry out the two subrac- 
tions (corresponding to mass and wave-function renormalisation) at q 2 = M 2 (“on-shell” 
scheme ||) i.e. 

dU(q 2 ) 


n R (q 2 ) = n (q 2 ) - n(M 2 ) - ( q 2 - m 2 ) 


dq 2 


(4.24) 


q 2 =M 2 


and so the renormalised self-energy II/j(g 2 ) becomes 


fl R (q 2 ) = ^rl (■f q 2 + ) ( L(q 2 ) - L(M 2 ) ) - ^( q 2 - M 2 ) ( 3 - 2 L(M 2 ) ) \ (4.25) 


47T 


1 

2 J 


11 

~8 


where 




(4.26) 


and A (q 2 ) was defined in 


Note that only the self-energy contribution II (q 2 


needs to be renormalised; indeed, after the PT rearrangement the resulting expressions for 


the vertices (and boxes) are ultra-violet finite, exactly as happens in normal QCD [23 


II Any other subtraction point p 2 would work equally well. 
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As a result, the gluon wave-function renormalisation constant and the gauge coupling 
renormalisation constant Z g are related by the QED-like relation Za = Z X J 2 |23|j2^ , |25| . 

In the limit q 2 M 2 , the leading (logarithmic) contribution to Il^/g 2 ) is given by 


n R (g 2 ) = 


2\ _ a s c A 9\ 2, / 2 / 2 


4tt + - 


(4.27) 


where the ellipsis denotes subleading contributions and /i is an arbitrary reference momen¬ 
tum. Instead, the corresponding limit for QCD is given by 


nrV) = ^ (y) tw//- 2 ) + ... 


(4.28) 


It is also interesting to compare the qualitative features of the MYM self-energy II 
with Cornwall’s massive propagator [Q which has been used successefully for fitting data 
T3| ; it has the form [^] (for Euclidean q 2 ) 

q 2 + 4 m 2 (q 2 ) 


d c {q ) = [q + m (q )}bg In 


A 2 


with 


2 / 2 \ 2 
m (q ) — m 


ln[ 


g 2 +4m 2 11 -12/11 
A 2 


, 4 m 2 . 


(4.29) 


(4.30) 


L ln (V 

where A is the QCD mass. Both d~^ l {q 2 ) and II (q 2 ) display the correct threshold behaviour 
(i.e. they turn imaginary for — q 2 = 4m 2 ). In addition (and in contrast to II(g 2 )) d^ 1 {q 2 ) 
has the correct asymptotic limit for q 2 A 2 , since the coefficient multiplying the leading 
logarithm is 11/3 (instead of 29/8 in the case of II(g 2 )) , thus capturing the one-loop 
QCD running coupling. Notice also the non-trivial dependence of the mass m(q 2 ) on the 
momentum. 

Finally, it is straightforward to check that if one inserts the expression for AmII(g 2 ) 
obtained from the tree-level calculation of the previous section into a twice-subtracted 
dispersion relation then one obtains the real part of the right hand side of (f4.25|), i.e. 


2\2 


$lefl R {s) = 


(s - M 2 ) 


7r 


4M 2 


ds' Amfl(s') 

(s' - M 2 ) 2 (s' - s) 


(4.31) 


We end this section by commenting on how the naive massive model gives precisely 
the same result for the one-loop S'-matrix element in question. We know that this must 


**The functional form for d c 1 (q 2 ) given in ( |f.29 ) represents an excellent, physically motivated fit to the 
numerical solution of a Schwinger-Dyson equation for the gauge-independent QCD gluon self-energy. 
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be the case, given the work of Section || The equivalence also follows from the tree-level 
arguments of the previous section. To see this, note that in computing the imaginary 
part of the one-loop amplitude we needed the qq —> gg amplitude only (i.e. in the unitarity 
equation the sum is over physical states and so ghost states do not contribute). In addition, 
the internal gluon propagators couple to conserved currents and so are the same in both 
MYM and naive gluon calculations. Under the assumption of analyticity, it follows that the 
two approaches give the same one-loop amplitude. To see how things go working explicitly 
with the full one-loop amplitude one needs to repeat the calculations of this section. The 
only differences between the S-matrix element of the MYM compared to the naive model are 
the replacement of the bare gluon propagator of ( |4.1|) with the unitary gauge propagator: 


Ufiv(q) I57ti/ 


%<lv \ 1 

M 2 ) q 2 — M 2 ’ 


(4.32) 


and the fact that the naive model does not have any ghosts. The actual calculation is 
straightforward, given the results presented above. One needs to replace the massless poles 
appearing in the auxiliary integrals Ji, J 2 , and Ji, stemming from the longitudinal part of 
the gluon propagator, by M 2 . The algebraic cancellations go through in exactly the same 
way as before with £ — ^ and - jp. 


5 Two-loop analysis 


Now we turn to the two-loop calculation. We will show that in this case renormalisabil- 
ity breaks down, and that the non-renormalisable terms are propagator-like and depend 
logarithmically on the cutoff. We will work again directly with the S-matrix element for 
the process qq —> qq. The calculations will be carried out using the Feynman rules for the 
naive massive gluon model since we know that, at the S'-matrix level, it is equivalent to 
the MYM. 


Consider the tree-level amplitude of Section |5| 7Y), (note that we have changed nota¬ 
tion by adding the subscript ‘0’ to denote a tree-level amplitude). It satisfies the following 
BRST identities [|6j: 


U^q-ab _ 

k2vS() b i 


l.uq-ab _ 

^2 2 0[iv ~ 

h. cab 


h^h u T ab 

M 2 S “ h , 

(5.1) 
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where 


(5.2) 


k a h a 

cab _ eabc n, l t re _ n rabc r ' , 2 t tc 

- 91 K ~ 9J V a . 

Using (|5.1|) one finds that the imaginary part of the amplitude, A, can be written: 

A = \j % ILU Q^(h)Q‘”’(k 2 )TS l „d(P.S?) 

= \ J {tTTL - S„Sl)d(P.S 2 ) , (5,3) 

whereas for normal Yang-Mills: 

Aqcd = t JP' v (h,v) p w °(h,n) T^d{P.S 2 ) 

= \ I - 2S„Sl,)d(P.S 2 ). (5.4) 

Notice that, despite the different factors accompanying the SoSq terms in ( [5.3|) and ( [5.4]) , 
both expressions give rise to renormalisable real parts, i.e. the real part can be ob¬ 
tained by means of a twice-subtracted dispersion relation. Renormalisability is manifest 
since the tree-level amplitudes contain no dangerous terms (such terms vanish by current 
conservation) and the contraction via the polarisation tensors does not induce any non- 
renormalisable terms, as a consequence of (|5.1|) . 

Proceeding to the two-loop analysis, one must consider two separate quantities: 
A2g(s, a 3 ), which is the contribution to the imaginary part of the two-loop amplitude 
which arises from the convolution of the tree-level amplitude for q{jpi)q{jp2) —> g{k\)g{k 2 ) 
with the hermitian conjugate of its one-loop partner, see Fig.il> and Az g (s, ct 3 ), which is the 
contribution to the imaginary part which arises on convoluting the tree-level amplitude for 
the process q{pi)q{j>2) g(ki)g(k 2 )g(k 3 ) with its hermitian conjugate, see Fig.[|. The two 

contributions must then be fed into a twice-subtracted dispersion relation and integrated 
from AM 2 to oo and from 9 M 2 to oo respectively, i.e. 


9£e fi(s, a 2 )) 


(s - M 2 ) 2 

71 


■AM 2 


ds' 9m Tl 2 p (s', a 2 ) 
(s' — M 2 ) 2 (s' — s) 


OO 



9 M 2 


ds' 3m n 3 ff (s', a 2 ) 
(s' — M 2 ) 2 (s' — s) 


(5.5) 


where 

A ng (s,a 3 ) = 2D 0 (s) Rrn fi ng (s,a 2 ). (5.6) 

If our calculations reveal that the RHS of (|5.5| ) is infinite then we will have shown that the 
5-matrix element for qq —>■ qq computed in the framework of the MYM is not renormalisable 
at two loops. 
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Figure 3: Diagrams contributing to 7^,. The blob refers to the same corrections to the 
gluon propagator as in Fig||a,e,f). 


We shall now show that -4. 2s , when fed into the integral on the RHS of (|5.5|) , gives a 
finite contribution, wheras the A^ g integral needs an additional subtraction in order to be 
rendered finite, i.e. it gives rise to a non-renormalisable contribution. 

To see that the contribution from *4. 2fl contains no dangerous terms, it suffices to 
prove that (a) the one-loop amplitude T° b u for qq —> gg is renormalisable, and (b) that 


i.e. 


it satisfies exactly the same type of BRST identity as its tree-level counterpart T^ b u , 
that (|01) holds if we replace T 0 ab u —> T^ u and —■> S^ b . Both (a) and (b) can be easily 

proved based on the analysis of [p8fl - It turns out that the closed expressions for T^ b v and 
Si b are given by the Feynman diagrams of regular QCD in the Feynman gauge, but with 
the tree-level propagators inside all graphs replaced by massive ones, again in the Feynman 
gauge, with the exception that for the ghost contributions we have a different statistical 
factor. This discrepancy does not affect the high energy behaviour of such graphs, i.e. the 
ghost loops are well-behaved for large q 2 . In addition, as follows from ^8j, the tree-level 
BRST identities do indeed hold at one loop. Thus 

1 


.4 


2 g 


(s,« 3 ) = - j2Me[T 1 ^Q^(k l )Q™(k 2 )T 0 l a ]d(P.S/ 
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Figure 4: Diagrams contributing to We do not show the diagrams which are related 
to those shown by exchange of outgoing bosons. 

= \ j 2 *e[(irisb - S,4)]d(J>.S. 2 ) (5.7) 

and the real part can be obtained using the twice-subtracted dispersion relation. 

Now we turn to the amplitude Tff p {k\, k- 2 , fc 3 ) for the process q{p\)q{p2) — > 

g(ki)g(k 2 )g(k 3 ), where fcj is the four-momentum of the z-th gluon, and pi + p 2 = q = 
k\ + k 2 + k 3 . Such an amplitude is given by the sum of the diagrams shown in Fig.[| 

Let us compute the quantity 

-^(s.a 3 ) = U K.W(fci) (5,8) 

where f d(P.S. 3 ) denotes the integration over the three-body phase space, with the com¬ 
binatorial factor 1/3! accounting for the three indistinguishable gluons in the final state. 
As before, the polarisation tensors satisfy the transversality condition: hq ■ Q(ki ) = 0. At 
first sight, the integrand in (|5.8|) seems to contain terms proportional to (M _2 )°, (M -2 ) 1 , 
(M -2 ) 2 and (M~ 2 ) 3 . The term proportional to (M -2 ) 0 is renormalisable, whereas all higher 
powers give rise to non-renormalisable contributions: The higher the power, the worse the 
divergence. However, as we shall shortly see, by virtue of the BRST identities that T pup sat¬ 
isfies and the transversality properties of the polarisation tensors, only terms proportional 
to (M -2 ) 0 and (M^ 2 ) 1 survive. Thus, the worst divergence is logarithmic. 

To establish this fact, let us first study the action of the longitudinal momenta k%, 
k%, and k 3 on T pup . It is straightforward to verify that Tff p {k\, k 2 , k 3 ) satisfies the following 
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identities: 


KT™ = (5 12 )f fe- + (S a n f , 

KT$Z = (Sn)f’k l „+{Sn)? 1 K, 

KT£, = {Su)t c h„ + (SsOffc, • 

Bose symmetry imposes the following relations among the S l3 amplitudes: 


(5.9) 


ss^fakj'kt) 

CLiCL£ 

— &ji 

Sij ajae (ki, kj, kf) 

^ di dgdj 

= 


and 


kj, ■ Sji — kj • S u , / i j■ 


The closed form of d>f| c (fci, k 2 , ks) reads 


\ abc 


(5 2 3 )T= (5|, 3 )r+(^) 


with 


r cs \abc _ „2 T /q f elc flab^l^aa^iki q, k\) f eal f lbc^3a(ki + k 2 ) 

\^23>u - 9 v e J J -?T7- 7 — T1777 -b / J 


elb rlac 


2k 2 ■ k 3 + M 2 


+ / 7 


(5.10) 

(5.11) 

(5.12) 


2a^3p 


(4j)f = -is 2 ®(P2) ( T“7„ ttTt r 7 „ + r 7 „ -A. 


2Ad • + Tf 2 2fci • k :i + M 2 

jeac j^a 


iX «(Pl) 


2fc 2 • /C3 + M 2 


(5.13) 


In deriving the above expressions, in addition to the elementary Ward identity, 
have employed the tree-level Ward identity: 


, we 


r S(?i> 92, 93, 94) = / a6 T$,(® > , 94, 9i + 92) + / ace r gf a (94, 92, 9i + 9s) 

+/ ade r^(92, 93, 9i + 94 ), (5-14) 

which relates the bare three- and four-gluon vertices. All remaining S l3 amplitudes can be 
obtained from S\ 2 using the relations of (|5.10|) . 

We next let the longitudinal momenta in the polarisation tensors act on T^ p and 
XjXt, and use (^9|). We can see how the (M -2 ) 3 terms disappear. The action of the term 
M^ 2 k^k^ gives terms proportional to k 3 and k 2 (or equivalently kf and k 2 ) , which vanish 


when they hit Q ,J/J {k\) or Q vX {k 2 ). So, (|5.8|) reduces to 


^3 g - — 


3! 




Xp 


d(P.S 3 ). 


(5.15) 
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At this point, the highest possible power of M 2 is (M 2 ) 2 . We now let M 2 /c 2 /c2 act on 
and Tj Ap : 

M- 2 k v 2 k$% vp Q^ih) rJ Xp = (s^t + sy^) Q pa (h) (s 2lp k° + 

= S&SgQ'vih). (5.16) 


Thus, one more power of M 2 has been eliminated. We are left with 

Aa, = ^j[(T m Tl t -S&&)<r’(,hMP-S.\ 

= Ar + Anr, 


(5.17) 


where 

Ag = 1 /FWtL - 3(S 12 Sj 2 )]d(P.S. 3 ), (5,18) 

A™ = I (*, ■ S 23 )(fc • S 23 )'d(P.S 3 ). (5.19) 

In deriving the above expressions we have used the fact that the phase-space integration is 
invariant under Q «->• kj. The term in (|5.19| ) will generate non-renormalisable terms; it is 
clearly non-vanishing, since it is a three-body phase-space integral over a positive definite 
quantity. In particular, 


a nr _ 

^g ~ 


r 2 

C A 


v:(h 


k 2 ) a (k 1 - k 2 fV*d(P.S 3 ) 


(5.20) 


3! 8M 2 

where the identities f abe f cde -f facefdbe _j_ jade jbce _ g an q jalmjbmnjcnl _ 1 CA j abc ] laye 

been used. We also used 2 kdc +M 2 = § + ■■■) where the omitted term is proportional to M 2 , 
thus giving rise to a renormalisable contribution, i.e. the omitted term belongs effectively 
to A^ g . Notice that the non-renormalisable terms are purely propagator-like (universal, 
process-independent). 


Finally, it is instructive to compare the result of ( [5.171 ) with that of normal QCD. In 
the QCD case there are, of course, no terms proportional to M~ 2 or higher powers. On the 
other hand, the presence of the auxiliary four-vector in the polarisation tensors could 
in principle induce spurious divergences, should it survive in the final answer. It is easy to 
see however how any reference to pg disappears before any of the phase-space integrations 
are carried out. Let us denote the corresponding QCD amplitude by T pU p. We start again 
with 


a qcd 

Aig 


^ /[Vp i) P fr (k 3 ) f} XT ]d(P.S 3 ) (5,21) 

where the gluons are now massless and /c* • P{ki ) =0 (7 = 1,2,3). Equation ( |5~7)| ) is valid 
for QCD, as can be shown rigorously using BRST arguments. Since the elementary Ward 
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identity ( |3.9|) is the same for both MYM and massless QCD, it follows that the closed 
expressions for the factors S tJ in Q5.9|) may be recovered from (|5.13D simply by setting 
M 2 = 0. We denote them by SV,-. Then, by letting the longitudinal momenta act on T pup 
and using (|5Tl), one can easily verify that any reference to the four-vector g p disappears 
and that the final answer is 

4f D W / - 6S a Sl 2 }d(P.S 3 ). (5.22) 

So, unlike the MYM model, in massless QCD all potentially dangerous terms vanish. 


6 Connection to field theories with a Higgs mecha¬ 
nism 


It is well known that the only way to endow gauge fields with mass whilst maintaining 
unitarity and renormalisability is via the Higgs mechanism [19,20,21 . This procedure 


is not suitable however for an effective model of strong interactions because it introduces 
extra scalar particles in the physical spectrum. In the MYM model massive gauge fields are 
obtained without introducing extra physical fields, at the price of losing renormalisability 
at higher orders of perturbation theory. It is instructive to see explicitly how the lack of 
renormalisability in the MYM can be traced back to the absence of a Higgs particle; in 
particular understand why the process qq —> qq (and qq —> gg ) is renormalisable at one- 
loop, but ceases to be renormalisable beyond one-loop. In this section we address these 
issues in detail by performing a quantitative analysis of the differences and similarities 
between the MYM and a Higgs model (HM) at the level of the 5-matrix. In addition, 


as has been discussed in detail in |L6|, it is possible to speak about the MYM using the 
language of a HM. Specifically, one can think of the MYM as a theory where all gluons 
have been given masses by adding to the Lagrangian a sufficient number of Higgs multiplets 
(N fundamental representations of SU(N ) ), and then “freezing” all the polar excitations 
of the N Higgs fields. The remaining N 2 — 1 angular excitations corresponding to the 
Goldstone bosons are precisely the angular fields 6 a displayed in (|2.1(J|). 


To illustrate the above points we will use a toy Higgs held theory which displays all 
the essential features we want to study. The gauge group of this model is SU(2). The Higgs 
mechanism is triggered by a complex doublet 0 in the fundamental representation (isospin 
l = |). This particular assignment endows all three gauge bosons with the same mass M, 
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whilst simultaneously prohibiting terms of the form 0-0-0 for any fermion representation of 
isospin i To mimic QCD, we choose the fundamental representation for the massless 
fermion fields -0, although this choice is not essential for what follows. As there are no 
interactions between fermions and scalars, the fermions remain massless even when the 
scalar fields acquire a non-vanishing vacuum expectation value. The Lagrangian density 
for this model is 

£ = - J V' 1 *' + + (D^)(Dyy - \/(0) (6.i) 

with 


(D,J>y = d^ + ig^jAl, 

1/(0) = /x 2 00 t + A(00 t ) 2 , (6.2) 


where T a = |cr a , and a a are the Pauli matrices. If fi 2 < 0, the Higgs mechanism gives 
rise to three degenerate massive gauge bosons of mass M — gv/2, where v = \J —pp is the 
minimum of V (0). The above model is a vector-like variant of the usual electroweak sector 
of the Standard Model, SU( 2 ) x 1/(1), with the Weinberg angle 6 W set to zero. 

The corresponding bare gauge boson propagator in the R% gauge has the form 


= 



-i 

q 2 — £M 2 )q 2 -M 2 


(6.3) 


and the would-be Goldstone boson (G) and ghost (c) propagators are 


A g, c = 


q 2 - £M 2 


(6.4) 


In the unitary gauge, which formally corresponds to the limit £ —> oo, the gauge boson 
propagator takes the form ( [4.32| ) and there are no Goldstone boson/ghost propagators. 0 
Finally, there is a Higgs particle of mass Mjj = v\/2X with bare propagator 


A H = 



(6.5) 


ttlt is elementary to verify that no gauge singlet (total l = 0) can be formed out of the above isospin 
assignments. 

ttThe renormalisability of the HM in the unitary gauge is not manifest. For example, it is known that, 
even though the n-point functions are non-renormalisable, by virtue of subtle cancellations, the S'-matrix 
element built out of these non-renormalisable n-point functions can be made finite with the usual mass 
and charge renormalisation [p2[ . A more immediate way to see this is to resort from the beginning of the 
calculation to the PT rearrangement of the amplitude [^9) . 
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Figure 5: Higgs boson contributions to . In the unitary gauge recall that (b) is absent. 


The one-loop /3-function for the gauge coupling has the form 


m 


p 


1 

167T 2 


{b g -b f - b s )g 3 


( 6 . 6 ) 


with 




b s 



(6.7) 


where Tj is the Dynkin index of the fermion representation, T s is the Dynkin index of the 
scalar representation, n/ is the number of fermion families in a given representation and n s 
the number of real scalar families. For the particular scalar representation we have chosen, 
T s = \ and n s = 2. In the absence of quarks we have that (3 = — jJ^-(43/6). 

Let us now proceed with a study of the one-loop amplitude. As we already showed, at 
the level of S'-matrix elements the MYM and the naive model are equivalent. In addition, 
if we adopt the unitary gauge for the HM then it is obvious that, to any finite order 
in perturbation theory, the only difference between an S'-matrix element computed in the 
MYM and the corresponding S'-matrix element computed in the HM is due to contributions 
to the latter which come from Feynman diagrams containing Higgs boson propagators 
f3l| . For example, in the case of one-loop quark scattering, in addition to the graphs 
in Fig.|2j. which are common to both MYM and HM, the diagram of Fig.[5|(a) contributes 
to the S'-matrix element of the HM. In other words, the S'-matrix elements of the MYM 
model may be obtained from the corresponding S'-matrix elements of HM by omitting all 
diagrams containing a Higgs particle. Given this observation, it is easy to see why the the 
process qq —> qq in the MYM is renormalisable at one loop: The only difference between 
the renormalisable HM and the MYM is the contribution corresponding to the graphs of 
Fig|| which themselves form a gauge-invariant and renormalisable subset. Denoting their 
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'WWW 


Figure 6: The Higgs graphs contributing to Tff. The graphs of Fig. [I] must also be included 
when computing T X H . 


contribution by Tiff we have that, up to the immaterial tadpole graph, 


fi»» = g 2 M 2 


d n k 

z(27t)’ 


~9ij, v + 


k^ky 


M 2 J (,k 2 - M 2 ) [{k + q) 2 - Mjj]' 


( 6 . 8 ) 


The factor M 2 in front of the integral originates from the gluon-gluon-Higgs coupling and 
guarantees that Tiff can be made ultraviolet-finite by means of the usual mass and wave 
function renormalisation. After (on-shell) renormalisation, the above expression becomes, 
in the limit q 2 M 2 , (and dropping the terms proportional to q^q u ) 

ft?" tf) = ~ (T) f In tf/M 2 )<T + ... (6.9) 

where the ellipsis denotes numerical constants and terms of order 0(M 2 /q 2 ). If we now 
set ca = 2 in the expression of (|4.27|) , and add it to the expression in (|6.9|) , we see that 
the coefficient in front of the resulting logarithmic term is equal to 43/6, which is precisely 
the coefficient of the HM /9-function without quarks. In addition, as expected from the 
discussion on the connection between the MYM and the HM given at the beginning of 
this section, the expression in ( |6.9| ) is exactly the difference between (4.27) and (4.28) for 
ca = 2. 


We now turn to the two-loop analysis. Firstly, it is relatively straightforward to estab¬ 
lish that the contribution from the two-gluon cut of those one-loop diagrams which contain 
a Higgs boson, Fig|| gives rise to a renormalisable contribution, i.e. the corresponding 
dispersive (real) part can be made finite by means of a twice-subtracted dispersion rela¬ 
tion. This is of course expected, since the one-loop diagrams contributing to A 2 g which we 
studied in the previous section were themselves renormalisable, i.e. no cancellation from 
diagrams containing a Higgs boson is required. 
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Figure 7: One of the graphs contributing to 7^ . The others are obtained by permuting 
the outgoing bosons. 


To see this explicitly, consider the amplitude A% g given by 


K„M = 5 / 2 SepljL Q"^) Q" a (k 2 ) T^ a }d(P.S.‘ 


( 6 . 10 ) 


where 7)^, is shown in Fig|| It is straightforward to verify that T^ w is a gauge-independent 
quantity, and that it satishes 

( 6 . 11 ) 


= h uS( 


H 


with 


s« = s WD"(,)[n;"B»( 9 )](fe - k 2 )“. (6.12) 

Notice that the expression in square brackets behaves like log(s/M 2 ) for s M 2 . Using 
fliTTl) and (| 6 . 11 |) we can see that 

1 


KM = 2 


2Se[(T 1 "T ( T t - S«s!,)]d(P.S 2 ). 


(6.13) 


So, this contribution gives rise to renormalisable two-loop amplitudes. 

To see how the presence of the Higgs boson enforces renormalisability, we focus on 
the two amplitudes, A^ g a and A^ g b : 

KfM = jJlT l S„Q'"(ki) QMMQ^iMKxMP-S 3 ) (6.14) 

and 

K" = | / 2Se[r"„ Q^(k,) KMP-S?) (6.15) 

where the amplitude T^ p is shown in Fig.[Fj. As can be seen, A^ g a arises by multiplying 
only those three-gluon amplitudes which contain Higgs particles whilst A^ g b comes from 
interfering the Higgs diagrams with the non-Higgs diagrams of Fig.^j. Since the coupling of 
the Higgs boson to two gauge bosons is proportional to M, it follows that 7^ has already 
a factor M 2 built into it. Consequently, there is an implicit factor M 4 inside A^ g °(s), and 
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therefore the only non-renormalisable contribution in M^’ a (s) will come from the term in 
the polarisation tensors which is proportional to (M -2 ) 3 . We therefore find that 


,11 1 

'2 3! M 2 


- fc 2 ) a (fc! - k 2 )%d{P.S 


(6.16) 


The interference term, A^g b , has an implicit M 2 inside so now contributions from the 
(M~ 2 ) 3 and (M -2 ) 2 terms in the polarisation tensors are needed. One finds, using (|5.9|) , 
the following non-renormalisable contribution: 

1 1 


[•A.3g b ]NR(s) = 


3! M 2 


v:(h - k 2 nh - k 2 fv 0 e d(p.s 3 ). 


(6.17) 


In arriving at the above results, identities of the type kfki = |(Aq — k 2 ) a (ki — k 2 ) 13 , or 
^ 2^3 = — ?^ 3^3 va hd under the integral sign, may be found useful. Comparing to (|5.20|) 
of the previous section (setting ca = 2), we see that the Higgs contribution exactly cancels 
the non-renormalisable part of the MYM two-loop contribution. Evidently, even though 
the Higgs boson does not couple directly to the quarks, (since in this toy model the gauge 
symmetry prohibits Yukawa couplings) its importance in restoring the renormalisability 
of the process qq —> qq manifests itself through the tree-level sub-amplitudes gg —> gg 
containing the Higgs boson (Fig.7), which reside in the two-loop diagrams. 


7 Quark-Quark Elastic Scattering 


In this section we take an introductory look at the elastic scattering of a pair of quarks 
via two-gluon exchange within the MYM model. Of course, quark-quark elastic scattering 
cannot be measured directly, but it is possible that many of the elements which are central 
to the more realistic processes (e.g. hadron-hadron elastic scattering) are contained in this 
simpler treatment. This is in the spirit of the Donnachie-Landshoff-Nachtmann approach 
■1-|3^ . where the success of the additive quark rule provides evidence that one need not 
know about the detailed structure of the colliding hadrons before one can proceed to make 
elastic scattering calculations, although it is not yet established that this is correct [[I3[ . 


As a first step, one can calculate the amplitude for the elastic scattering of differently 
flavoured quarks, i.e. q t q 3 —> q t q 3 ^ at the lowest order, keeping only those terms which 
dominate in the Regge limit. This is a straightforward calculation of the box diagram 
shown in Fig.|(c) rotated through 90 degrees (the crossed box diagram contributes only to 
the real part of the amplitude in the Regge limit and constitutes a sub-leading correction). 
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The leading contribution is imaginary and so can be obtained directly using the cutting 
rules, i.e. the amplitude for single gluon exchange can be written: 


A(s, = iT <8 T°2ri k2 9 ^ p 2rt6 Xl 


(7.18) 


where k is the momentum of the exchanged gluon and pi and p 2 are the momenta of the 
incoming quarks, i.e. s = (p\ + p 2 ) 2 . The high-energy limit allows the exchanged gluon 
to be assumed soft, and so the eikonal approximation has been used to simplify the qqg 
vertex. The delta functions ensure hclicity conservation at each vertex. Multiplying by the 
conjugate amplitude, projecting out the colour singlet part and performing the two-body 
phase space integral (putting the intermediate quarks on-shell) allows us to write 


A(s,t) 


isa 2 


N 2 - 1 
N 2 


d 2 k 


1 

(k 2 + M 2 )((k-q ) 2 + M 2 )’ 


(7.19) 


where q 2 = — t > 0, N is the number of colours and the exchanged gluons are taken to be 
purely transverse. The transverse momentum integral can be performed and yields 


s _ 2 W 


„ .... 11 Alt) +1 

A(s, t) = i —27ra;—— -—— In ——--, 

v ’ ' -t s N 2 A(t) A (t) -1’ 


(7.20) 


where A(£) is defined in (|3.28|) . Thus, the total cross-section for q t q 3 —> X is 


,N 2 - 1 


ar ~ Jp™* 


N 2 


(7.21) 


It is instructive to investigate the conditions under which the two-gluon exchange 
amplitude calculated above violates unitarity. We shall see that unitarity is violated only 
for very central collisions and that these constitute an insignificant fraction of the total 
and elastic scattering cross-sections. Only for very high-f processes do we have collisions 
which are sufficiently central to cause a worry. This gives us confidence to proceed to the 
next order of calculation, assured that we have yet to receive indications that unitarisation 
corrections are important. 

To investigate unitarity we perform a Fourier transform of the elastic scattering am¬ 
plitude, i.e. 

~ sb )~ f d2q c — 

J( 2ir) 2 2s ’ 

and b is the impact parameter of the collision. Written in this way, unitarity demands that 

|lM)| 2 <i 
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for all b. However, we can be confident that nnitarisation corrections are small if the 
inequality is satisfied for those values of impact parameter which dominate the process 


under study. Numerical evaluation of (|7.22|) demonstrates that the amplitude only ever 
violates unitarity for Mb < 10~ 2 , 10 -3 for 2/a 2 = 20, 50 respectively, i.e. only for very 
central collisions (on the scale of the gluon mass). In this language, the total cross-section 
is given by 

a T = 2 J d 2 b A(s, b) (7.23) 

whilst the elastic scattering cross-section is given by 


(Tei = / d 2 b \A(s,b)\ 2 . 


(7.24) 


Since A(s, b) decreases monotonically as Mb increases, it follows that the elastic scattering 
cross-section receives a larger contribution from more central collisions than the total cross- 
section. To a first approximation, the typical impact parameter is set by the gluon mass, 
i.e. (b 2 ) = C/M 2 where C r-u 1 and is larger for the total cross-section than for the elastic 
cross-section. In either case, we are always well away from the dangerous region where 
unitarity is violated. The situation is different for high-7 processes, since now (b 2 ) 1/1*1 

and so, for yf—t > (10 2 — 10 3 )M we would need to worry that unitarisation corrections are 
important. 


8 Conclusion and Perspectives 

In this paper we have reviewed and investigated the formalism of the MYM model, ar¬ 
guing that it may be relevant as a tool to investigate diffractive scattering (and possibly 
other areas of strong interactions phenomenology), where traditional QCD methods are 
inadequate. A detailed study of the qq —> qq process in the context of this model up to 
the two-loop order was presented, and the renormalisation properties of the corresponding 
S'-matrix were discussed. 

Let us summarise briefly the prospects for a study of diffractive scattering. A suc¬ 
cessful model of diffraction should be able to explain: The growth of total hadronic cross- 
sections with increasing s. In particular, the model should show why the rise in soft 
hadronic processes (e.g. the total pp cross-section) proceeds at a much slower rate than in 
hard processes (e.g. the cross-section); the shrinkage of the forward diffraction peak 
with increasing s. In other words, the model should be able to explain the qualitative 
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success of the Donnachie-Landshoff-Nachtmann model of the pomeron as a single Regge 
pole in soft diffraction (i.e. in those processes where there is no large scale) and its failure 
in small-x deep inelastic scattering, in the diffractive production of all vector mesons at 
high Q 2 and in the diffractive production of J/'h mesons at low Q 2 . 

In the future, we plan to use the MYM model to compute the complete 0(a s ) cor¬ 
rections to the two-gluon-exchange amplitude discussed in the previous section, in order to 
verify whether some or all of the aforementioned features emerge. More specifically, such 
a calculation should help us investigate the following points: 


In the limit of large enough s, the logarithms ~ (a s lns) m , which appear at each 
order in perturbation theory, become large and it becomes necessary to sum them to 
all orders. This summation of leading logarithms is performed using the formalism 
of BFKL. It is an open question as to precisely when this summation leads to the 
dominant contribution to the amplitude. In fact, since the summation is of leading 
logarithms only, we cannot define exactly what it means to say s is large, although 
we note that, in this respect, analysis of the next-to-leading logarithmic corrections 


calculated by Fadin, Lipatov et al. |34[ should improve the situation. By computing 
at fixed order in a s we can investigate the relative importance of the ~ a s In s term 
compared to the terms which do not include the logarithm. In this way, we can make 
some quantitative statements regarding the need (within the MYM model) to sum 
the remaining leading logarithms. For example, it might be that, at the energies of 
contemporary colliders, the logarithm is not so large to justify dropping the other 
terms, i.e. a fixed order calculation might be the better way to proceed. 


It is known that introducing a gluon mass has a very small effect on the leading 
logarithmic contribution |35].|4[]. This arises largely because the BFKL summation is 
infrared finite, i.e. there are infrared cancellations between real and virtual graphs 
which reduce the sensitivity to this region. These cancellations persist even after 
adding a gluon mass (via a Higgs, or via the MYM model) and serve to reduce the 
sensitivity of the amplitudes to variations in the mass. Note that is not too important 
how the mass is introduced. This can be seen since no Higgs graphs contribute to 
the leading logarithm summation and since the gauge-dependent part of the gluon 
propagator is also sub-leading (in covariant gauges). 


• It is also known that the leading logarithm summation leads to a rapid rise of total 
cross-sections. It can be argued that this rapid rise, which is due to multiple soft gluon 
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emission, reveals itself in hard scattering processes, but is masked in softer processes 
by unitarity corrections. Any slowing down of the rise via unitarity corrections has 
yet to be precisely quantified. Another possibility is that the strong rise seen in hard 
processes can be explained in fixed-order perturbation theory, i.e. arising from the 
In s term, and that this same rise is masked in soft processes by a non-logarithmic 
contribution which is comparable in size to the logarithmic contribution (i.e. as the 
process becomes harder, the non-logarithmic contribution falls away to reveal the 
logarithm). This latter possibility can be investigated after computing the radiative 
corrections to the two-gluon exchange graphs. 

In summary, we think that the MYM may prove a useful tool in understanding the phe¬ 
nomenology of diffractive scattering by bridging the gap between different QCD-inspired 
models. Such a conjecture will be tested through next-to-leading order calculations of 
quark-quark elastic scattering, which we plan to discuss in a forthcoming paper. 
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